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Abstract 

The systematic approach for the calculations of the non-perturbative contributions to the free energy 
in the ferromagnetic phase of the random field Ising model is developed. It is demonstrated that such 
contributions appear due to localized in space instanton-like excitations which exist only in dimensions 
D < 3. It is shown that away from the critical region such instanton solutions are described by the set 
of the mean-field saddle-point equations for the replica vector order parameter, and these equations can 
be formally reduced to the only saddle-point equation of the pure system in dimensions (_D — 2). In the 
marginal case, D — 3, the corresponding non-analytic contribution is computed explicitly. Nature of the 
phase transition in the three-dimensional random field Ising model is discussed. 
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1 Introduction 

Despite extensive theoretical and experimental efforts during the past several decades very little is understood 
about the basic thermodynamic properties of the ferromagnetic Ising systems with quenched random fields 
(for reviews see e.g. [1]). In the most simple form the random field spin Ising systems can be described in 
terms of the Hamiltonian 

N 

F = - ^ CTjO-j - ^ hiffi (1) 

<i,j> i 

where the Ising spins {ai = ±1} are placed in the vertices of a D-dimensional lattice with the ferromagnetic 
interactions between the nearest neighbors, and quenched random fields {hi} are described by the symmetric 
Gaussian distribution: 



P[h^] = (27rhl) exp 



1 



N 

^ l—l 



(2) 



where the parameter hg describes the strength of the random field 

According to simple physical arguments by Imry and Ma [2] it should be expected that the dimensions 
Dc above which the ferromagnetic ground state is stable at low temperatures (it is called the lower critical 
dimension) must he Dc = 2 (unlike the pure Ising systems where Dc = 1) Indeed, if we try to test the stability 
of the ferromagnetic state by flipping the sign of the magnetisation in a large region of linear size L, we will 
find two competing effects: a possible gain of energy, due go alignment with the random magnetic field, which 
scales as Eh = hoL^^^; and the loss of energy, due to the creation of an interface, which scales as L^~^. 
These estimates show that below dimensions Dc = 2 for any non-zero value of the field at sufficiently 
large sizes L the two energies are getting comparable, and therefore no spontaneous magnetization should be 
present. On the other hand, at dimensions greater than Dc — 2, the energy of the interface is always bigger 
than Eh- Therefore these excitations will not destroy the long range order, and a ferromagnetic transition 
should be present. These naive (but physically correct) arguments was later confirmed by a rigorous proof 
by Imbrie [3] . 

On the other hand, a perturbative study of the phase transition shows that, as far as the leading large 
scale divergences are concerned, the strange phenomenon of a dimensional reduction is present, such that 
the critical exponents of the system in the dimension D appear to be the same as those of the ferromagnetic 
system without random fields in the dimension d=D-2 [4]. Since the lower critical dimension of the pure 
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Ising model is equal to 1, this result would imply that the lower critical dimension of the random field Ising 
model must be equal to 3, in contradiction with the rigorous results. Actually, the procedure of summation 
of the leading large scale divergences could give the correct result only if the Hamiltonian in the presence of 
random fields has only one minimum. In this case the dimensional reduction can be rigorously shown to be 
exact, by the use of supersymmetric arguments [5]. 

However, one can easily see that as soon as the temperature is close enough to the putative critical point 
(as well as in the whole low temperature region), there are local values of the magnetic fields for which the free 
energy has more than one minimum [6] . In this situation there is no reason to believe that the supersymmetric 
approach should give the correct results and therefore the dimensional reduction is not grounded. Thus, the 
above arguments settle the controversy about the lower critical dimension of the random field Ising model in 
favour of the value Dc = 2. 

Although at present we understand that at D > 2 the low temperature state of the RFIM must be 
ferromagnetic, the nature of the phase transition from the paramagnetic to the ferromagnetic state is still 
a mystery. The only well established fact here is that the upper critical dimensionality (the dimensionality 
above which the critical phenomena are described by the mean-field critical exponents) of such systems is 
equal to 6 (unlike pure systems, where it is equal to 4). What is going on in the close vicinity of the phase 
transition, at dimensions D < 6 is not known. And it is not just the question of what are the true values of 
the critical exponents. Even the nature of the phase transition (is it of the second or of the first order) is nor 
clear. Moreover, at least in some cases there are indications in favour of the existence of the intermediate 
spin-glass state (with broken replica symmetry) separating paramagnetic and ferromagnetic phases [7, 8] . 

The key problem here is the effect produced by the non-perturbative contributions which arc coming from 
numerous local minima states. In the present study, considering these contributions in the low temperature 
ferromagnetic phase (away from the critical point) at the finite strength of random fields ho, I would like 
to "rehabilitate" the marginal character of the dimensionality D = S as well as of the trick of dimensional 
reduction D {D — 2). The point is that although numerous local minima configurations always exist 
in the ferromagnetic phase of the RFIM, they are not always statistically relevant for the thermodynamic 
properties of the macroscopic system. In this paper I am going to demonstrate, first in terms of simple 
heuristic arguments (Section 2), and then using more rigorous replica technique (Section 3), that the non- 
perturbative states yield (non-analytic) contribution to the thermodynamics only in dimensions D < 3. In 
terms of the qualitative heuristic arguments this contribution appears due to rare large cluster spin flips. In 
terms of the non-perturbative replica formalism [9, 10] such configurations are described by the instanton type 
solutions of the mean-field saddle-point equations for the replica vector order parameter with broken replica 
symmetry. These equations can be formally reduced to the only saddle-point equation of the corresponding 
pure system in dimensions (£> — 2) (this fact has been noted first in [11]), and then one can easily note that 
localized instanton- like solutions formally exist only in dimensions D < 3. In this sense the dimension D = 3 
is marginal: in dimensions D > 3 the non-perturbative states are irrelevant, while in dimensions D < 3 
they yield finite non-analytic contribution (finally, when approaching the dimension D = 2 from above these 
states become so much relevant that they destroy the ferromagnetic ground state of the system). Similar 
instanton-like configurations in the presence of the homogeneous external magnetic field have been recently 
studied in Ref.[12]. Although formally the saddle-point equation in dimension £) = 3 (in the zero external 
field) have no localized instanton solutions, nevertheless, as it often happens in the marginal situations, the 
the relevant contributions can be taken into account in terms of the instanton-like configurations containing 
the soft-mode parameter (which is the size of the instanton). 

Away from the critical point the concentration of the instantons (or the flipping clusters) is exponentially 
small, so that thoy can be treated as independent finite energy excitations. However, when approaching 
the putative phase transition point (from below) the typical distance between instantons eventually becomes 
comparable with their size. In this situation the present scheme of calculations breaks down, as the state of 
the system becomes just a mess of locally ordered "up" and "down" regions, when it is impossible to separate 
the degrees of freedom connected with the flipped clusters from those of the ferromagnetic background. It is 
shown however, that at the finite strength of the random fields this happens at temperatures where the system 
can still be described at the mean-field level, and it is argued that in this situation it would be reasonable to 
expect the phase transition of the first order into the disordered state (Section 4). 
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2 Heuristic arguments 



2.1 Perturbative contributions 



To compare the perturbative and the non-perturbative effects in the random field Ising model let us consider 
its Ginzburg-Landau continuous representation: 



if[0(x) 



\ (V0(x))^ + ir02(x) + - /i(x)</.(x) 



(3) 



Here r = {T/Tc — 1) (|r| <C 1) is the reduced temperature parameter (for simplicity in what follows it will be 
supposed that = 1) and the random fields h{x) are described by the Gaussian distribution, 



P[h{y 



■ po exp 



is/"" 



X /i2(x) 



(4) 



where ho is the parameter which describes the strength of the random field, and po is an irrelevant nor- 
malization constant. The system will be considered in the low-temperature ferromagnetic phase, so that the 
reduced temperature parameter r will be taken to be negative, t = — \t\. To neglect the effects of the thermal 
fluctuations (away from the critical point), it will be assumed that the absolute value of the parameter |t| is 
not too small: 



Irl » o^/^^-^) = 



TGL 



(5) 



which is the usual Ginzburg-Landau condition. The only relevant spatial scale in the system described by 
the Hamiltonian (3) is the correlation length: 



Rcir) 



-1/2 



(6) 



(which under condition (5) is described by the mean- field critical exponent v = 1/2). The statistical prop- 
erties of the system at scales bigger than the correlation length can be studied in terms of the saddle-point 
configurations defined by the equation 

-A0(x) - |t|<^(x) + 503(x) = /i(x) (7) 

In the absence of the random fields the ferromagnetic ground state of the system is described by the 
homogeneous configuration (j)o = ^/\t\/ g, and it has the free energy density /o = — T^/(4(7). In the usual 
perturbative approach the effects produced by the random fields can be taken into account e.g. in the following 
way. After rescaling the fields 



0(x) = ( ^ 



1/2 



instead of eq.(3) one obtains the system which is described by the rescaled Hamiltonian 



H 



0(z) 



(V0(z)) 



. (i-p) 



V9 



where z = x/Rc = |t|^/^x. Here rescaled random fields 

h{x/Re) = R-^ 
are described by the distribution function 



|x'-x|<flc 



d^x' h{x') 



P[h{z)] =poexp(-^^^2^^ J d^'zh^z) 



(8) 



(9) 



(10) 



(11) 



which is characterized by the mean square value 
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= \Tf''hl (12) 
The ground state configurations of tlie fields ^(z) are now defined by the saddle-point equation 

-A^(z) - 0(z) + 4>\2) = ^h{^) (13) 

In the absence of the random fields {h = 0) the above equation would yield the trivial homogeneous ferro- 
magnetic ground state solution ^(z) = 1 (or ^(z) = — 1). In the presence of the random field term the ground 
state solution can be represented as 

^g«(z) = l + ^(z) (14) 

where V'(z) describes small {ijj <C 1) spatial fluctuations of the ground state configuration. Solving eq.(13) 
in the perturbative way (considering its r.h.s term as the small perturbation), for the typical value of these 
spatial fluctuations one easily finds: 

r-,A + o((,^)'] (15) 



|t|(6--D)/2 ^ 1 \^|t-|(6-D)/2 

Substituting the solution eq.(14) into the Hamiltonian eq.(9) for the free energy density of the ground state 
configuration one gets: 



1 r~ 1 \r\^ I ~\ 1 IrP 



4g V 'JR^ Ag 



where 5f is the random quantity with the typical value 



9K , ( 9K 



+ o ) ) (17) 



We see that the perturbation expansion goes in powers of the (small) parameter 

which means that such perturbative approach is bounded by the condition 

\r\:$>T.{h,)^{ghl)^ (19) 
On the other hand, according to eq.(9), the typical value of the thermal fluctuations 

(20) 

In the pure system wc can neglect the effects of the thermal fluctuations provided {dP") ^ 1, which imposes 
the usual Ginzburg-Landau condition, eq.(5). Hero, in the presence of random fields, the thermal fluctuations 
can be neglected provided they are small not only compared to the average value of the ferromagnetic order 
parameter (which in the present notations is equal to one), but also compared to the typical value of its 
random spatial variations (which is described by the fleld ij), eq.(14)): 

(^')«V^ (21) 
Substituting here eqs.(15) and (20) we find one more bound for the temperature parameter |t|: 

|r| « hi (22) 

Thus, taking into account eq.(19), we conclude that the perturbative mean-field consideration of the present 
system is legitimate in the temperature interval 

2 

,2\ 6^ 



ghiy-- ^\T\^hl (23) 
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which automatically requires that the typical value of the random field can not be too small: 



hi > (24) 

Note that under this condition 

TGL-^T^iho) (25) 

Another words, when approaching the phase transition from below (which means reducing the absolute value 
of the temperature parameter |t|), first one arrives to the crossover temperature t*(/io) when the random 
fields can not be considered as small perturbations any more (while the thermal fluctuations can still be 
neglected), and only after that, at |t| ~ tgl (where the perturbation theory is no more valid) the thermal 
fluctuations would become important. 

In the particular case of the dimension £> = 3 (which is the main focus of our further study) the expansion 
parameter of the perturbation theory is 



|r|3/2 



« 1 (26) 



The crossover temperatures are 



TGL = 



r{ho) = {ghlf (27) 
and the consideration is limited by the temperature interval 

T{ho) «: |t| « hi (28) 
which can exist provided the typical value of the random fields is not too small: 

/lo > (29) 

Within these bounds the perturbative calculation of the free energy can be expressed in the form of eqs.(16)- 
(17). 

2.2 Non-perturbative contributions 

Besides perturbative contributions described above, one can observe completely different type of thermal 
excitations. Let us suppose that in the low temperature phase the system has the global ferromagnetic 
magnetization directed "up", and let us consider a spatial island with sufficiently large linear size L, in which 
the average value of the (random) field h is negative and sufficiently strong. Then one can easily note that 
in addition to the state "up" (with slightly modified value of the order parameter), another local minimum 
with orientation "down" can be formed within this island. Considering our system in terms of the rcscalcd 
Hamiltonian, eq.(9), one can easily see that this alternative state can exist only if the gain in the energy due 
to the interaction with the random field, 

Eh ~ (30) 
overruns the loss of energy due to the creation of the interface, 

Ef ^ +\IL^L^-^ (31) 

g 

These estimates demonstrate that for a given value of L such double-state situation in the considered island 
is created provided 

\h\ > K{L) ^ LJ_ (32) 
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which has a finite probability for any finite value of L. 

Let us consider the statistical properties of such type of excitations in more detail. In terms of the 
Hamiltonian, eq.(9), the energy of the spherical region of the radius L in which the magnetization is flipped 
"down" can be represented as the sum of two contributions: 

E{L) = (const) ' L°-^-V{L) (33) 
.9 

The first term here is the positive energy due to creation of the spherical interface, while the second contribu- 
tion is the random energy which appear due to the interaction with the random field, which in the continuous 
limit it can be represented as 

V{L) = 2 '^' ' / d^zh{7) (34) 

V5 i|z|<L 

According to this definition, 

I |(l-£)) I I <^-°> 

V^{L) ~ L^h? ~ ^ hi (35) 



Since the typical value of the random energy 



1/2 



scales as i^/^ we conclude that in dimensions D > 2 

the first term in the total energy E{L), eq.(33), must dominate at large scales, and therefore E{L) is (on 

average) a growing function of L (this is nothing else but the familiar Imry-Ma arguments [2] which explain 
why the ferromagnetic state is stable at D > 2). The point, however, is that E{L) is the random function, 
and it grows with L only on average, while for a given realization of disorder it can have local minima at 
various (large) values of L. 

The probability distribution which describes the statistics of the random functions V{L) is given by 



V[ViL)] = Jvhiz) P[hiz)] n J^^^^^d^zh{z)-V{L) 



(36) 



where P[h{z)] is the probability distribution function, eq.(ll). Performing straightforward Gaussian integra- 
tions one obtains 



■P[V{L)] = (const) exp 



dL 



(dV{L)y 



(37) 



where Sd = 2t^"I^ /T{D /2) is the square of the unite sphere in D dimension. Since we are going to consider 
only largo scales, the above result is given in the continuous limit representation containing the ultraviolet 
cut-off of the order of the correlation length, which in the present notations is equal to one. We see, that in 
accordance with the physical meaning of the function V{L), its statistical distribution depends only on its 
derivative (the constant term in this function is irrelevant), and therefore the problem would become well 
defined only if we fix a value of this function at a given L. For simplicity, let us assume that at the lowest 
possible scale V(L = 1) = 0. 

The above probability distribution function, eq.(37), can be used to estimate the probability that the 
random function E{L) has at least one local minimum at scales larger than a given scale L. Since the value 
of the energy E{L) in a putative minimum growth with L, its probability must be small at large scales. In 
this situation the sufficient condition for the existence of a minimum somewhere beyond a given size L is 
^<0,or 

^ > (con.t)^^L-- (38) 
dL g 

The probability Pmin{L) that the above condition is satisfied at a unit length at the given size L can be easily 
estimated using the general distribution function (37). Formally it can be obtained by integrating P[F(L)] 
over all functions V{L) conditioned by eq.(38). It is clear, however, that with the exponential accuracy, the 
result of such integration is defined by the lower bound for the derivative dV{L)/dL. Therefore, substituting 
(at a given value of L) ^^^f^ = (const) i|T|' 2 '1,^-2 Jj^^q eq.(37), with the exponential accuracy one gets: 
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Pmin{L) ~ exp 



9% 



(39) 



Thus we have derived very important property of the random function E{L), eq.(33): although, according 
to cq.(35) this function (in dimensions _D > 2) on average grows with L, the probability of finding a local 
minimum of this function in dimensions D < 3, according to eq.(39), also grows with L. The situation in 
the three-dimensional case is marginal, and it is not quite clear to what extent the above simple arguments 
are grounded for _D = 3 (as usual in the marginal situations, more rigorous methods has to be used). 
Nevertheless, if we formally substitute D = 3 into eq.(39), we would have to conclude that the probability of 
a local minimum becomes independent of the size of the flipped cluster. Since the value of Pmin {L) \ ( d=3) is 
exponentially small (in the parameter \t\'^^'^ / {gho) ^ 1, eq.(26)), the contribution to the free energy of such 
type of cluster excitations with the exponential accuracy is defined by their probabilities (while their energies 
define a pre-exponential factor). Thus we can estimate the non-perturbative part of the free energy density in 
the ferromagnetic phase of the 3D random-field Ising model as 

13/2 



A/ ~ exp -(const) ^-^-2- (40) 
V 9% J 

In the next section we shall re-derive this result in terms of the systematic replica approach, which, in 
particular, allows to calculate the {const) factor. 



3 Non-perturbative replica approach 

The general scheme of the non-perturbative replica calculations has already been discussed in detail in the 
recent papers [9, 10]. Here we repeat it just in brief. Let us consider a general D-dimensional random system 
described by a Hamiltonian H [4){x)-. h{x.)], where (^(x) is a field which defines the microscopic state of the 
system, and /i(x) are quenched random parameters. Let us suppose that in addition to the ground state, 
this system has another thermodynamically relevant local minima states located "far away" from the ground 
state and separated from it by a finite barrier of the free energy. In other words, it is supposed that the 
partition function (of a given sample) can be represented in the form of two separate contributions: 



(41) 



where Fq is the contribution coming from the vicinity of the ground state, and Fi is the contribution of the 
local minima states. Then, for the averaged over disorder total free energy we find: 

1^ ^ 1: 



J^=--lnZ = Fo--ln[l + Z,Z^'] (42) 

The second term in the above equation, which will be denoted by AF, can be represented in the form of the 
series: 

1 (_T\m—l 1 °° (^'\Y 

\ ^ \ '-I rrm'T-m _ 1- \ ^ V ^) 



AF = -- V ^ Zf^o"™ = -- hm V ^ Z„(m) (43) 

where 



Zn{m) = 11 [ Vcj>^^^ n / e-''^4^^^'--^^^'-<''-'^n-".] (44) 

6=1 c=l '' 

is the replica partition function (if„ [0] is the corresponding replica Hamiltonian), in which the replica 
symmetry in the n-component vector field 0a (a = 1, n) is assumed to be broken. Namely, it is supposed 
that the saddle-point equations 

^=0, ,..V.„„, (45) 
have non-trivial solutions with the RSB structure 



7 



m = 



ii(x) for a = 1, m 

60 (x) for a = m + 1, . 



(46) 



with 01 (x) ^ 00 (x), so that the integration in the above partition function, eq.(44), goes over fluctuations 
in the vicinity of these components. It should be stressed that to be thermodynamically relevant, the RSB 
saddle-point solution, eq.(46) should be localized in space and characterized by a finite space size R.{m,) and 
finite energies E{m) = Hn [</>*]. In this case the partition function, eq.(44), will be proportional to the entropy 
factor V/R^{m) (where V is the volume of the system), and the corresponding free energy contribution AF, 
eq.(43), will be extensive quantity. 

Thus, to compute the non-perturbative free energy contribution one should find all saddle-point RSB 
solutions 0a (x), eq.(45) and their corresponding energies E{m) (in the limit n — > 0), and finally one has to 
sum up the series 



AF 



V ^ 



-R-^im) ((3deiT] 6 

V J n=0 



m 



Here 



5^H [0] 



<i>=<t>' 



(47) 



(48) 



and the Hessian prc-cxponential factor (dctT)^^/^ (taken in the limit n — > 0) appears due to the integration 
over the replica fluctuations in the vicinity of the RSB solutions, eq.(46). Note that in the present approach 
the procedure of analytic continuation n — > is quite similar to that in the usual replica theory [14]: whenever 
the parameter n becomes an algebraic factor (and not the summation parameter, or the matrix size, etc.), it 
can safely be set to zero right away. 



Coming back to the original Hamiltonian, eq.(3), and following the standard procedure, after the Gaussian 
averaging of the replicated partition function over the random fields /i(x), one obtains the replica Hamiltonian 



TL Ti n n 

L a=l a=l a=l o,6=l 

The saddle-point configurations of the fields (f)a{x) are defined by the equations 

-A0„(x) - |t|0„(x) -h50^(x) - hi Ij^M^)] = 



(49) 



(50) 



^5=1 



Below it will be demonstrated that besides the obvious (replica symmetric) ferromagnetic solution 0a (x) = 
4'o = \fV\f9 these equations have non-trivial localized in space instanton-like solutions with the RSB struc- 
ture: 



a:(x) 



^Vi(x-v/M) for a = l,...,m 



(51) 



^ Vo(x-\/|t|) for a = m + 1, ...,n 



These solutions arc characterised by two non-trivial functions V-'i(z) 7^ "02(2) (where z = x-\/|t|). Substituting 
these rescaled fields into the saddle-point eqs(50) and into the Hamiltonian, eq.(49), we find that (in the limit 
n — > 0) the instanton configuration {'0i(z), '0o(z)} is defined by the equations 



-A^i - ^1 V? - A(m) (V'l - ^0) = 

-A^o - ^0 + V'o - •^M (V'l - ^0) = (52) 

and its energy is 
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I |2-D/2 r 

E{m) = m^-^ J d'^z 



(53) 



where 



A(m) = ^ (54) 

We are looking for the localized in space (spherically symmetric) solutions of the eqs.(52), such that the two 
functions ^i(r) and tpo{r) (where r = |z|) are different from each other in a finite region of space, and at 
large distances they both sufficiently quickly approach the same value V'o = 1) so that the integral in eq.(53) 
will be converging. 

Keeping in mind the qualitative arguments of the previous Section it will be assumed that the parameters 
of the model satisfy the requirements, eqs.(23)-(25). In this case the considered theory, eqs.(52)-(54), contains 
the large parameter A(m) ^ 1 (for any m = 1,2, ...), so that, according to eqs.(52), the two fields Vi and V'o 
must be close to each other. Redefining, 

V'o(r-) = ij{r) - jx{r) (55) 
in the leading order in <C 1 instead of eqs.(52) we get much more simple equations: 

-Alp -tp + tp^ -2x = 

-Ax+{3^^-l)x = (56) 

which contain no parameters. For the energy of the configurations described by the two fields ■)/'('') and 
x{r) instead of eq.(53) (again, in the leading order in A~^) we find the value, which does not depend on the 
summation parameter m, 

where 

Eo = J d'^z [(VV')(Vx) + (V' - V')X - X'] (58) 

is the universal quantity which depends only on the dimensionality of the system. Considering only spherically 
symmetric configurations, eqs.(56) can be represented as 

D-ld^p „ 
^+V^-2x = 



dr"^ r dr 
d-'x D-ldx , 



dr'^ r dr 

Now as a matter of a simple algebraic exercise one can easily check that taking 



{3^P^-1)X = (59) 



X = --f (60) 
r dr 



the above two equations (59) can be reduced to one equation 



9 



for the only function ^p{r). The corresponding energy, eq.(58), of the configurations described by eqs.(60) 
and (61) can be reduced to 



(62) 



Thus, we have reduced the problem of the non-perturbative excitations in the D-dimensional random field 
Ising model to the study of similar instanton-like saddle-point configurations in the corresponding pure system 
in dimensions (D — 2). 

Physically relevant solutions of the saddle-point equation (61) must be such that in the limit r ^ oo 
the function ip{r) sufficiently quickly approaches the value Vo = 1 (so that the integral in eq.(62) would be 
converging). It is evident, that at D > 3 (which corresponds to the "effective" dimension Dp/f = (D — 2) > 1) 
eq.(61) has no such solutions. The formal solutions of this type and their contributions to the free energy in 
dimensions 2 < D < 3 (which correspond to < D^ff < 1) has been first considered in Ref.[ll], and then 
studied in detail in Ref.[10]. 

Here we are going to concentrate on the marginal case I? = 3, when eq.(61) reduces to 

^ = _ ^ (63) 
(where < r < +oo). Although this equation has no instanton solutions, the configurations of the type 



V(r) = tanh (^-^J (64) 

(which are the formal solutions of eq.(63) for — oo < r < +oo) at large enough values of the parameter L 
can be considered as the soft- mode "quasi-instantons". Using eq.(62) one can easily find that at L :» 1 the 
energy of such configuration is weakly dependent on the soft-mode parameter L: 

Eo{L) ~ - 8^277 exp(-4L/V2) (65) 

where 

-B* = —^'^ (66) 

is the energy of the infinite-size configuration. 

It is clear that the energy Eo{L) monotonously decreases with L, so that the true minimum is achieved 
when the above quasi-instanton configuration completely disappears. Therefore the configurations described 
by eq.(64) represents continuous spectrum of thermal excitations described by one parameter L. Correspond- 
ingly their contribution to the free energy is given by the summation (the integration) over all possible sizes 
L. In terms of the original spatial notations (x = z|r|~^/^ = zRc) the size of the above quasi-instanton is 

R = LRc (67) 

Thus, coming back to the general expression for the off-perturbative part of the free energy (with D = 3 and 
/? = 1) eq.(47), (where the energy E(rn) and the size R are defined in eqs.(57) and (67)), and introducing 
here the integration over the soft-mode parameter L, we obtain 

^ ^ ^ — ' ' ' -Eo{L) (68) 



AF = -i fcii^T^E^^expfJ 



ho9 



(here the ratio V/ (LRc)^ is the entropy factor which yields the number of positions of the instanton with the 

linear size [LRc) in the three-dimensional volume V). Note that for the derivation of AF with the exponential 
accuracy, in the considered range of parameters, eqs.(23)-(25), the contribution of the fluctuations (contained 

^ 1/2 

in the term (detT)^^g , eq.(48)) can be neglected [10]. Substituting here eq.(65), and neglecting all pre- 
exponential contributions, for the density of the non-perturbative part of the free energy we finally get (cf. 
eq(40)) 
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|r|3/2 



A/ ~ exp(-£;.^) (69) 



where = ^^^tt. Note that the vaHdity of this result is Hmited by the condition: (^/iq)^^^ ^ kl ^ '^o- 



4 Discussion 

The present study of the non-perturbative phenomena in the random field Ising model has been done in terms 

of the mean-field local minima configurations in the low-tempcraturc ferromagnetic state. The approach 
completely neglects thermal fluctuations, which is justified, on one hand, provided the temperature is not 
too close to the putative phase transition point, and on the other hand, provided the temperature is not too 
low. where (weak) thermal fluctuations would nevertheless overrun small spatial quenched fluctuations due to 
the interaction with random fields. In terms of the Ginsburg-Landau Hamiltonian, eq.(3), in the dimension 
D = 3 these two requirements impose the restrictions on the value of the reduced temperature parameter, 

^ |t| ^ ft-Q, which automatically implies that the considered procedure of calculations is formally valid 
only at finite strength of random fields, ho ^ g. Under these restrictions both the " hand- waving" heuristic 
arguments (Section 2.2), and the formal replica calculations (Section 3) provide the results, eqs.(40) and (69), 
which nicely fit each other. Of course, the main thing here is not the actual value of the non-perturbative part 
of the free energy, but the physical mechanism, which provides it. According to the speculations of Section 
2.2 the non-perturbative contributions appears due to rare large cluster spin flips (which in terms of replica 
calculations are described by the localised in space instanton conflgurations). It is crucial that such clusters (or 
the replica instantons) are supposed to be far from each other, so that they can be treated as non-interacting 
and independent. According to the obtained results, eqs.(40) and (69), the spatial density of these clusters, 
p ^ cxp[—{const)\T\^^^ / (ghg)], is exponentially small provided |t| '3> {gh^Y^^ = T{hQ). In other words, when 
approaching the putative phase transition point from below (i.e. decreasing the value of |r|), at |t| ~ T{ho), 
the mean separation between clusters becomes comparable with their typical size. In this situation the whole 
scheme of calculations breaks down, as the state of the system becomes just a mess of locally ordered "up" 
and "down" regions, when it is impossible to separate the degrees of freedom connected with the flipped 
clusters from those of the ferromagnetic background. It is crucial that this happens at the temperature 
t(/io) 3> TQL = .9^ (well before the putative Tc of the expected ferromagnetic/paramagnetic phase transition) 
which is still far way from the Ginzburg-Landau crossover temperature tgl where the thermal fluctuations 
would become important, and therefore here the system can still be described at the mean-field level. Since 
at |t| ^ t(/io) the local order parameter (the absolute value of the local magnetisation) is still finite, it would 
be reasonable to expect that somewhere at these temperatures the system undergoes the phase transition of 
the first order into the disordered state. As for the nature of this disordered state, it should be stressed that 
(unless the replica symmetry is broken!) one should not expect one more phase transition from the spin-glass 
to the paramagnetic state. The point is that due to the presence of quenched random fields the quantity {(j))^ 
(which is the traditional replica-symmetric spin-glass order parameter) remains non-zero at all temperatures. 
Therefore one should expect not more than a crossover from " rather spin-glass" (just after the transition) to 
"rather paramagnetic" (at high temperatures) disordered states. 

Unfortunately the analytic technique developed in this paper can not be directly applied for the description 
of this (rather exotic) disorder induced first-order phase transition. Nevertheless, since the situation seems 
to stay at the mean-field level, development of another analytic approach to this problem does not look 
completely hopeless. 
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